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Abstract We compare the results obtained from analytical perturbation theory and 
the AKM numerical code for an axistationary spacetime built from matching a ro- 
tating perfect fluid interior with the equation of state £ — 3p = 4B of the simple 
MIT bag model and an asymptotically flat exterior. We discuss the behaviour of the 
error in the metric components of the analytical approximation going to higher or- 
ders. Additionally, we check and comment the errors in multipole moments, central 
pressure and some other physical properties of the spacetime. 



1 Introduction 

The lack of stellar models in General Relativity -i. e. a stationary and axisymmetric 
perfect fluid interior matched with an asymptotically flat vacuum exterior- is in 
direct contrast with the importance they could have for the astrophysics of compact 
stars and in particular the determination of their possible compositions. One of such 
possibilities is the interesting case of strange matter. A versatile model for strange 
matter is a simple MIT bag model with equation of state (EOS) £ — 3p — 4B. We will 
study this stellar model with the results provided by the CMMR post-Minkowskian 
and slow rotation approximation scheme [3] and its behaviour when we go to higher 
orders of approximation. Also, we will give the relative error in these functions and 
quantities when compared with very precise numerical results obtained with the 
AKM code [1,2]. 
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2 The analytical and numerical metrics 

The stationary and axisymmetric spacetime we study is the following. The interior 
y~ is filled with a perfect fluid with the EOS £ — 3p = £q and in ridig rotation 
so that, if % are the Killing vectors, its velocity is u = \\f(% + with \j/ a 
normalization factor and CO a constant correspondig to the rotation speed of the 
fluid as seen by a distant observer. The exterior Y + is asymptotically flat vacuum 
and is matched with the interior imposing continuity of the metrics and their first 
derivatives on the p = surface. 

In CMMR, we solve the Einstein equations using a truncated multipolar post- 
Minkowskian approximation in spherical-like coordinates associated to harmonic 
ones. The post-Minkowskian parameter is A =m/r s , whith m the Newtonian mass 
of the source and r s the coordinate radius of the static fluid; a different parameter 
Q 2 = co 2 r 3 /m (the ratio between Newtonian centrifugal and gravitational forces), 
gives the truncation point of the expansion in spherical harmonics, in this case pre- 
serving terms up to &(Q 3 ). The metric in each spacetime g ± (A, Q) is decomposed 
in Minkowski TJ plus the deviation A ± (A, Q) and then Einstein's equations can be 
solved iteratively. The spacetimes are matched on the p = surface which can be ex- 
panded as r% = r s [l + oQ 2 i^fcos 0)] + &(Q 4 ) with a constant. Finally, the global 
metric depends only on £q, and r s . In [4] (CGMR) we obtained the <^(A 5 / 2 , Q 3 ) 
metric for the EOS e + (1 — k)p = £Q. Here we use its (k = 4, £q = 4B) subcase corre- 
sponding to the simple MIT bag model but now including terms up to <^(A 9 / 2 , Q 3 ). 

AKM is a multi-domain spectral code that gives the matched metric on a grid over 
a quadrant of finite size. The grid coordinates are {p, £} , cylindrical associated to 
quasi-isotropic coordinates and the resolution is customizable. It also gives some 
physical properties, like the first mass and angular momentum multipole moments 
Mo and J\, baryon rest mass Mb, circumpherential radius R c hc, binding energy E\>, 
polar and equatorial coordinate radii r p , r e , polar redshift z v and central pressure 
and specific enthalpy p Cl h c . To build a stellar model, it needs goal values for any 
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Fig. 1 Relative error for constant density in g tt using Mo = 8 x 10~ 4 , CO = 0.2 using CMMR up to: 
(a) <^(A 5 / 2 , Q 3 ); (b) ^(A 9 / 2 , Q 3 ). The thin dotted lines represent the AKM and CMMR surfaces 
(indistinguishable in this picture size) 
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Fig. 2 Relative error between ^(A 9 / 2 , Q 3 ) CMMR and n = 20 AKM in (a-b) g tt and g t(p for 
M = 8 x 10" 3 , 0) = 0.24 (A w 0.077, Q « 0.059); (c-d) g„ and for M = 0.0184, co = 0.24 
(A « 0.12, i2 « 0.059). The thin dotted lines represent the AKM and CMMR surfaces. 

pair of these quantities as well as to fix some parameters to specify the EOS and 
provide an initial data file. The code is able to obtain initial data for different EOS if 
one manages to avoid unphysical configurations in the process. The precision it gets 
depends on the number n of Chebyshev polynomials used, and can reach machine 
accuracy for high enough n when the deformation of the source is not extreme. We 
fix it here to n = 20. 



3 Comparison and results 

We build stellar models for different values of (Mb, co) and compare the metrics on 
the cylindrical-like coordinate grid of AKM. Working in units where (c = G = B = 
1), the k = 4 CGMR has only (r 5 , co) as free parameters. Unlike (0, r s has no equiv- 
alent in AKM, so we must adjust its value. Hence, for each model we take the value 
of one of the AKM physical quantities and equating it to its CMMR counterpart we 
get r s . Using Mo, p Cl and J\ for this gives similar results. Here we use Mq. 

The results of Fig. 2 show the relative error in the metric (g tt , g t § ) components for 
Mo = 0.5 1 M and the typical Mo = 1 .39M for a moderately fast rotation frequency 
V = 102 Hz near the source. The first case have a reasonable error (higher in g t ^) 
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and the second must be improved. Going to higher orders in the A approximation 
is now automatic with the help of our Mathematica subroutines, but improving the 
slow-rotation approximation is more cumbersome. Nevertheless, Fig. 1 shows the 
improvement for a £ = £o EOS of moving from ^(A 5 / 2 , Q 3 ) to ^(A 9 / 2 , Q 3 ). It 
make us confident that, since our current error graphics do not show the lobular 
aspect of Fig. 1(b), which means that the truncation of the multipolar expansion is 
an important error source, our results can be improved very easily even in the case 
of the strong gravitational field of a compact source of realistic mass. 



Table 1 Some CMMR values and relative error with respect to AKM in (c = G = B = 1) units. In 
other units, the three models rotate at v = 102 Hz and their masses in M appear in the first row. 



M 


M = 0.054, M b = 0.06 


M = 


0.51, M b =0.6 


M = 


1.39, M b = 1.7 




CMMR 


error 


CMMR 


error 


CMMR 


error 


(0 


0.24 




0.24 




0.24 




M b 


7.9997e-4 


3.2e-5 


7.95e-3 


6.1e-3 


0.02241 


0.0406 


Mo 


7.104e-4 




6.772e-3 




0.0184 




J\ 


8.3609e-8 


le-4 


3.586e-6 


0.012 


1.82e-5 


0.034 


^circ 


0.034769 


3.6e-5 


0.07287 


8.6e-3 


0.102 


0.068 


^bind 


8.8730e-6 


2.8e-4 


3.766e-4 


0.040 


1.75e-3 


0.19 


Zp 


0.021133 


7.6e-5 


0.1079 


0.016 


0.24 


0.13 


V 


-0.020913 


4.6e-5 


-0.103 


0.010 


-0.221 


0.087 


r v 


0.033914 


4.7e-5 


0.0657 


9.8e-3 


0.0831 


0.09 


r e 


0.034055 


3.8e-5 


0.06593 


9.8e-3 


0.0834 


0.091 


f ratio 


0.995862 


9.4e-6 


0.9964 


2.4e-5 


0.997 


3e-4 


£ c 


4.5935 


6.8e-5 


4.712 


0.019 


5.56 


0.21 


Pc 


0.043679 


2.4-4 


0.254 


0.052 


0.631 


0.37 


h c 


0.908786 


1.7e-6 


0.9522 


2e-3 


1.025 


0.042 
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